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We investigate the effect of non-degenerate quarks on Bk- This effect is noticeably large for Bk 
(significantly larger than statistical uncertainty). Hence, it is important to include this effect in 
order to determine Bk with higher precision. We also observe that the quality of fitting for Bk 
gets better when we include non-degenerate combinations to fit to the prediction by Van de Water 
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1. Effect of non-degenerate quarks on Bk 

This paper is a follow-up from the previous paper [jl]| on Bk- Hence, in order to avoid repeating 
the same explanation, we will quite often refer to Ref. [|l|] for details. 

In this paper, we focus on Bk, the kaon bag parameter for indirect CP violation. In particular, 
we are interested in the effect of non-degenerate quarks {i.e. ^ m,,). Before we go into the 
details, we want to address couple of issues which covers the historical evolution of the staggered 
e' /e project. When we used unimproved staggered fermions to calculate weak matrix elements, we 
have observed three major problems: (1) large scaling violation, (2) large perturbative corrections 
at the one loop level, and (3) large uncertainty due to the quenched approximation. We use im- 
proved staggered fermions with HYP fat links which reduce the scahng violations and decrease the 
one-loop correction down to ^ 10% level [§,|3[ §]. We measure weak matrix elements over a subset 
of the MILC gauge configurations in order to remove the uncertainty originating from quenched 
approximation [||]. In other words, we numerically simulate the unquenched QCD where the va- 
lence quarks are HYP staggered fermions and the sea quarks are AsqTad staggered fermions (we 
call this a "mixed action"). Hence, in this paper, we report recent progress in calculating Bk in 
unquenched QCD using improved staggered fermions. 

Details of input parameters for this numerical study are provided in Ref. [|l]] and so we do not 
repeat them here. We use 5 degenerate quark combinations and 4 non-degenerate combinations to 
probe the effect of non-degenerate quarks. These 9 combinations are listed in Table |l|. Note that 
there are two pairs of which the degenerate and non-degenerate combinations share the same kaon 
mass at the leading order (for example, = niy = 0.03 and = 0.01 &my = 0.05). 





Ifly 


Bk 


0.01 


0.01 


0.6108 ±0.0164 


0.02 


0.02 


0.6421 ±0.0065 


0.03 


0.03 


0.6722 ± 0.0039 


0.04 


0.04 


0.6955 ± 0.0029 


0.05 


0.05 


0.7144 ± 0.0024 


0.01 


0.05 


0.6869 ± 0.0068 


0.02 


0.05 


0.6907 ± 0.0039 


0.03 


0.05 


0.6976 ± 0.0030 


0.04 


0.05 


0.7058 ± 0.0027 



Table 1: Quark mass combinations and corresponding Bk values 

In Fig. [l], we present Bk as a function of Euclidean time with a degenerate quark combination 
(nix = my = 0.03) in the lefthand side and with a non-degenerate quark combination {nix = 0.01 
and My = 0.05) in the righthand side. The best fitting range is 10 < f < 15 as explained in Ref. ||l|]. 
Here, note that the quark masses are chosen such that the kaon masses are the same at the leading 
order for both the degenerate and non-degenerate combinations. In Fig. we observe two things: 

(1) the signal for the non-degenerate combination is noisier than that for the degenerate one and 

(2) the Bk value of the non-degenerate combination is noticeably higher than that of the degenerate 
one. 
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Figure 1: Bk vs. T for degenerate valence quarks with m^ = niy = 0.03 (left) and for non-degenerate valence 
quarks with nix = 0.01 and w, = 0.05 (right) 



In Ref. [||], Van de Water and Sharpe calculated the chiral behavior of Bk using staggered 
chiral perturbation. They provide results for the case of Nf = 2 + 1 partially quenched QCD (m„ = 
ntci ^ nis) in the continuum as well as at a finite lattice spacing with staggered fermions. For the 
degenerate combination (nix = my), the chiral behavior of Bk is given in Ref. ^ and we do not 
repeat it here. For the non-degenerate combination {rux 7^ niy), they predict the chiral behavior of 
Bk in the continuum limit as follows: 

/If 2 {^\~^yY' 2 2\ 

Bk = c\y\ + ^^^2p |_^™™ + ^disc + C2mxy + C3 h £4 {2mjj + m^j 

+ -A (1-1) 



where / = 132 MeV and c,- are unknown dimensionless low-energy constants. The chiral loga- 
rithms are given by /conn and /disc- The contribution from the quark connected diagrams is 

2 2 
/conn = (>m%l{m%) - ?>l{ml) + 3limj) (l + ^) (1.2) 



The contribution from the diagrams involving a hairpin vertex is 

F+/r,)/m^ 



l&sc = {h+lY+h)/mly (1.3) 



7 2 



-l{ml] 
-l{ml] 



{m\. + ml) {ml - ml) {ml - m|) 
(m2 -m|)2 

2{m\, + m\){ml, - m\){ml - m| 
(m2 -mj)(m2 -m|) 
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2. {ml- m|)(m| - m|) - {ml + m|)(m| - m|) - (m^. + m|)(m2 - m|) 



mi, — mt 



Iy = lx{m\ ^mj) 



Irj = l{m^ 



^^{ml-m\f{m\. + m\){ml-ml) {mj - m\ ] 



{m\ — m^)^{mY — m^Y 

Here, note that the chiral logarithms appear through 1{X) and 1{X), which are defined as 

1{X) = XIog(X/A2)+F.V. 
1{X) = -[Iog(X/A2) + l]+F.V. 



(1.4) 
(1.5) 
(1.6) 



(1.7) 
(1.8) 



where F.V. represents a finite volume correction. The notations for the various meson masses are 
as follows for those composed of sea quarks: 



(1.9) 



mj;=2/xmrf, m^ = 2/lm.v, m^ = (mu -\-2m^)/3 

and as follows for those composed of valence quarks 

m|- =2lJ.mx, my =2}Xmy, m^y = }x{mx + my) (1-10) 

Here, we set sea quark masses to m„ = 7^ and the two valence quark masses are m^ and my. 
0.8 



0.4 



"T 1 1 \ 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 — 

(mixed, tree, non— degenerate) 




J I I I I I I I I I I I I L_ 



0.1 



0.2 



0.3 



0.4 



0.5 



M/ (GeV^) 



Figure 2: Bk vs. 

In Fig. ^ we plot Bk data as a function of M^. Here, it includes both degenerate combinations 
(wijc = m^,) and non-degenerate combinations {mx / my). We fit the data to the form of Eq. (|0]) 
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suggested by chiral perturbation theory: 

1 



ci 1 + 



48712/2 



'conn ~l~ ^disc 



2 (m|-my)2 . 
+ C2m^, + C3 2 Vc^m^, 



m 



(1.11) 



xv 



where the cut-off scale A is set to A = Anf. The fitting results are summarized into Table |^ 



parameters 


average 


error 


c\ 


0.4469 


0.0159 


Cl 


-1.4629 


0.1788 


C3 


0.0216 


0.0075 


C4 


1.4206 


0.2133 


X^/d.o.f. 


0.2999 


0.2668 



Table 2: Fitting results for Bk 

When we compare the fitting results of Table ^ with those in Ref . (only for the degenerate 
combinations), we notice three things: (1) is much smaller for the whole fitting (including 
both degenerate and non-degenerate combinations) than that only for the degenerate combinations, 
(2) the low-energy constants ci, C2 and C4 are consistent with those of Ref. [|l]] within statistical 
uncertainty, and (3) for a given choice of A, the C3 term is negligibly small. Judging from these 
observations, we understand that the effect of non-degenerate quarks are dominated by the chiral 
logarithmic terms such as /conn and /disc whereas the C3 term does not play much role at all. Judging 
from the fact that the C3 term is so small, we know that we can determine all the low-energy 
constants (ci, C2, C4) using only degenerate combinations. In addition, using these low-energy 
constants, we can determine the Bk value at the physical kaon mass with physical valence quarks 
(nix = rrici and ruy = rus). However, we know that is substantially smaller when we include non- 
degenerate combinations in fitting, which leads us to the conclusion that it would be significantly 
better for data analysis to measure enough non-degenerate combinations as well as degenerate 
combinations. 

In Table ^ the x'^ is rather high when we take into account the fact that all the data are 
correlated. In fact, one can observe that the fitting curve does not fit the lightest two data points 
very well in Fig. In other words, the fitting curve miss the lightest two data points in the opposite 
direction. In Ref. [^, it is pointed out that the contribution from the non-Goldstone pions are so 
significant that the curvature of the fitting curve becomes smoother, which is consistent with what 
we observe in Fig. ^ However, the full prediction from staggered chiral perturbation contains 21 
unknown low-energy constants for a single lattice spacing and 37 unknown low-energy constants 
for the full analysis [|6|]. In order to determine all of them, it is necessary to carry out a significantly 
more extensive numerical work. 

In the current analysis of Bk data, we match the lattice results to the continuum values at the 
tree level. In this respect, the results are preliminary. 

2. Effect of non-degenerate quarks on 5^^^^^ and 5g^^^^ 



In Fig. |[, we plot B^^^'^'' in the lefthand side and b'^^''^^ in the righthand side as a function of 
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Figure 3: fi^^''^^ vs. M| (left) and fif vs. M|. (right) 



M|^. We show data of both degenerate and non-degenerate quark combinations. The difference in 
Bf'^^ (fif between the quark mass pairs of (0.03, 0.03) and (0.01, 0.05) is about 0.6% (0.7%), 
which is much smaller than that in Bk (about 2%). Hence, the effect of non-degenerate quarks on 

(3/2) 

By^ is of less interest as much. 

There ai^e some calculation on the chiral behavior of {n'^ \Q^^J^}^ available in Ref. |^. This 
result applies to the Nf = 3 case with rriu = = rris. It, however, does not apply to the case of 
Nf = 2 + 1 with niu =md ^ m^. Therefore, in order to understand the data in Fig. |3| further, it is 
necessary to calculate the chiral behavior for Nf = 2 + 1 using staggered chiral perturbation theory. 
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